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Abstract 

We present a one-loop calculation of a gauge invariant QCD beta function. Us- 
ing both momentum and temperature renormalization group equations we investi- 
gate the running coupling in the magnetic sector as a function of temperature and 
momentum scale. At fixed momentum scale we find that, in contrast to A</)^ or 
QED, high-temperature QCD is strongly coupled, even after renormalization group 
improvement. However, if the momentum scale is changed simultaneously with tem- 
perature in a specified manner, the coupling decreases. We also point out in what 
regime dimensional reduction occurs. Both the cases Nf smaller and larger than 
^^Nc are discussed. 
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1 Introduction 

One of the most difficult aspects of finite-temperature QCD, and in fact finite-temperature 
field theory in general, is the appearance of large perturbative finite-temperature correc- 
tions at temperatures high compared to any zero-temperature mass scale. This may be 
illustrated for the case of theory in 4D. There one finds for the two-point function 



where A is the zero-temperature coupling and m is the zero-temperature mass. Clearly 
the above expansion is ill-defined in the regime T ^ m. One possible way out is to 
perform a partial resummation of the above as a diagrammatic series, i.e. summing the 
ring or plasmon diagrams. This is basically equivalent to expanding around a theory of 
mass ~ AT^ rather than a theory of mass m^. If this is done one finds a non-analytic 
expansion in A ||T| 



In critical phenomena language the above is akin to looking at a mean field theory 
or large- mass limit, the large mass being an induced finite-temperature mass mx. Such 
partial resummations would be inadequate to investigate the vicinity of a phase transition, 
such as the deconfinement phase transition or the electroweak phase transition. Near a 
second-order or weakly ffist-order phase transition one expects to be expanding around a 
very small effective-mass theory as opposed to a large effective-mass theory. Because in 
this case T 3> mx we can regard this regime as being a high-temperature regime as well. 

These problems encountered in finite-temperature field theory arise because one tries 
to describe one asymptotic regime in terms of effective degrees of freedom associated with 
another. If one wishes to describe finite-temperature systems, one needs to take into 
account that the effective degrees of freedom in the system are scale and temperature 
dependent. Generically there will be (at least) three asymptotic regimes wherein the 
physics is qualitatively different. Again these can be simply illustrated in a A0^ theory. 
The ffist one is when T/m ^ oo with m the zero-temperature mass. In this case one gets 
driven into a large (finite-temperature) mass regime. This mass subsequently acts as an 
IR cutoff thereby eliminating the IR problem. The second is when T/mx — > 0, mx — > 0. 
This corresponds to the massless 4D theory. The third is when T/mx — > oo and this 
corresponds to the massless 3D theory. 

There are fixed points of a temperature-dependent renormalization group (RG) associ- 
ated with all three asymptotic regimes as shown in [Q. This corresponds to the fact that 
in the three different regimes the effective degrees of freedom that describe the theory 
are qualitatively different being infinitely massive 4D, massless 4D or massless 3D re- 
spectively. If one implements an appropriate temperature dependent RG then the entire 
crossover between high- and low-temperature behaviour becomes accessible. The results 
show that A(/)^ exhibits a second-order phase transition and that the critical exponents 
associated with this phase transition are 3D ones, being associated with a 3D fixed point, 
not a 4D one. 

In this note we wish to implement this philosophy to a further understanding of finite 
temperature QCD. In section 2 we present the ingredients of a one-loop calculation, and 
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derive the expression for the /3 function. In the hterature there has been ambiguity 
to even the sign of the (3 function due to gauge, vertex and configuration dependence. 
Landsman has shown how to get around this by using the background field method 
in the Landau gauge which is a special case of the Vilkovisky-de Witt construction of a 
gauge-independent effective action [^]. In section 3 the momentum and finite-temperature 
(5 functions are integrated to give the finite-temperature running coupling. Finally in 
section 4 we discuss our results. 



2 Renormalization 

Because Lorentz invariance is broken at finite temperature, the QCD self energy in the 
Landau gauge decomposes into three- dimensionally transverse and longitudinal parts. 
Being independent, the transverse and longitudinal parts may behave differently under the 
RG, and also different masses ("electric" and "magnetic" respectively) may be generated. 
In the static limit (zero external energy) for a heat bath at rest, the (inverse) propagator 
only has non-vanishing time-time and space-space components. The transverse part 11*'' 
of the self energy is then found from the trace over the spatial components 



n-dPl) = ,3) 



Owing to gauge invariance, the background field method yields additional Ward iden- 
tities that give rise to relations between the scalar structure functions in the tensor decom- 
position. In particular the spatial, static three-point vertex in the symmetric momentum 
configuration can be connected to the transverse self energy and the tree vertex in the 
following way |^ 

rkVt^M I ^ = 9 r'^n [g,,{p, - q,) + cycl.] [l + n*=-(|p |)] . (4) 

The diagrammatics are specified by the special background-field QCD Feynman rules of 
Rebhan . These rules are easily extended to finite temperature 0, 0] . 

It is well-known that renormalization of the theory at zero temperature is sufficient 
to make Green-functions ultraviolet-finite at all temperatures. However, we have already 
mentioned that at high temperatures, if one uses perturbation theory in terms of zero 
temperature parameters, one gets large corrections that make the finite-order perturba- 
tive results unreliable. Therefore one expects a better description of the temperature 
dependence if one expresses the theory at finite temperature T in terms of parameters 
characteristic of that temperature. The essence of this is that when one sets a renormal- 
ization condition, one is deciding to parameterize the physics using the parameter defined 
by the condition in question. 

We renormalize the theory by requiring as a renormalization condition that the spatial, 
static, three-gluon vertex function is equal to the tree vertex at the momentum scale k, 
and temperature r. Formula (^ then implies that 

n*^(b-'U,T,r)||^^l^,^^^, = 0. (5) 

Note that there are now two arbitrary renormalization scales k, and r. The assumption of 
renormalization point independence of the bare vertex functions gives rise to a momentum- 
scale RG and the finite temperature RG [§]. Using eq. (||) one finds that the corresponding 
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beta functions for the effective coupling as{K, r) = r)/47r^ are 



dn 
d 



d\p\ 
dli^' 



r^a,(/t,r)=/?W(/t/2r) = «, T 



\p\=K,T=r 



(6) 
(7) 



\p\=K,T=T 

where a factor of two in the argument of the (3 function is included for later convenience. 



3 Results 

In the one-loop self-energy calculation we use the retarded/advanced formalism developed 
by Aurenche and Beccherawy [^], which gives directly the analytic continuation from the 
imaginary Matsubara energies and should in the static case give the same result as the zero 
external Matsubara frequency. Because of the presence of distribution functions in the 
diagrammatic rules, conventional methods like dimensional regularization and the use of 
Feynman parameters are less convenient than in the zero-temperature case. For instance 
in a one-loop integral it is not always possible, by shifting the integration variable, to 
remove angular variables simultaneously from the argument of the distribution function 
and from the propagators. And if we use Feynman parameters, when we change the 
Feynman parameter from zero to one, the pole of the integrand can move from below 
to above the real axis in the complex energy plane. We have chosen to regularize the 
one-loop self-energy integral in four dimensions by letting the differentiation with respect 
to the external momentum act on the integrand of the loop integral and by keeping the 



ie in the retarded/advanced propagators finite [ p!0[] . 

After summing the diagrams and performing the angular integrations the final result 
for the beta functions can be expressed in terms of the following one-dimensional integrals: 



F:i{y) 
Gl{y) 



dx x"log ^^-^ — 

e^y — r] X — 1 

dx P — — — , 

e^y — 1] [x^ — 1)" 



(9) 



where P stands for the principal part. For non- negative n the integrals converge at both 
ends. The result is: 



F 

vao 



The bose thermal part and vacuum part are given by {y = k,/2t): 

3^-,_ 21_-,_ 25_i_ 7 
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B 

vac 



6 



(10) 
(11) 



(12) 
(13) 
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The fermion parts are found to be (if we neglect fermion masses compared to temperature 
and momentum scale): 



F 

vac 



Nfal 
\Nfal 



f-,--F^\y)--Fo\y) + Gi\y) 

QV 4 4 



(14) 
(15) 



Expansions of the integrands around zero show that these thermal expressions are infrared 
convergent. The distribution functions take care of the ultraviolet convergence of the 
thermal parts and also make these parts vanish in the low-temperature limit. 

Derivatives of the functions and can be expressed in terms of these function 
themselves by performing an integration by parts 

l^y^ + 2m+ljFUy) = 2 G?(y) + 2 J] — [l - (1 - r^)^-^'^] C(2fc) r(2fc), (16) 



y^G'liy) = 2Gl{y). 
dy 



(17) 



Using these identities we can write the thermal parts of the beta function as a derivative 
of a function / 



+ i^l = ^lyi-f{y) 



d 



dy' 



18) 



where 



f{y) = 

+Nf 



IT 



21 ..1... , 7^1, 



12,^ ^-,^^iy)^V6^^iy)^-sGliy) 



1 „ , _ 1 



2Ay^ 4 



+ F,\y) + -F,\y) 
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(19) 



This function / vanishes in the low-temperature limit. 

The fermion parts are in complete agreement with calculations in the literature 0] 
(with conventional and with background-field Feynman rules; the one-loop fermion con- 
tribution is the same in both cases). In the boson term, however, there are small numerical 
differences with the results of Antikainen et. al. |^ which may be due to the fact that 
they used a background field method in the Feynman gauge. 

Integration of the beta functions from a point {kq, tq) in the k, r-plane with value 
asi^o, To) as initial condition, then gives the running behaviour of the coupling as a func- 
tion of the momentum scale and the arbitrary temperature r so that it makes the vertex 
function invariant of the renormalization scale up to a wavefunction renormalization. Thus 
the "experimental" input here, as{i^o,To), is the coupling at the temperature tq and at 
some fiducial momentum scale kq. From this one describes all the physics at other scales 
and other temperatures. 

Integration of the system (^ and (0) of coupled differential equations gives: 
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-N^ 



In h 



/ 



K 

27 



-1 



(20) 
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By rewriting the integration constant asi^o, ^o) as 
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Nf]\n 



A(/to, To) 



(21) 



where now A(ko, Tq) becomes the arbitrary parameter, we can bring the above expression 
for the running couphng ( PO] ) in the form 
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Nr - -Nf In 



A(/€o,ro) \2t 



[22) 



which is useful for comparison with vacuum calculations of the running coupling constant. 
In the zero-temperature limit the function / vanishes, so that we recover the standard 
vacuum QCD one-loop result with A(ko, 0) = Aqcd the usual vacuum QCD scale. As 
the renormalization scale r is arbitrary, we can freely choose r equal to the physical 
temperature T and in the next section we will do so. Similarly the momentum scale k is 
set equal to the physical momentum \p\. 



4 Asymptotic behaviour 



For ?/ — > oo we find for the (3 functions 



(23) 



which is the standard /3 function for QCD in 4D with massless fermions, and Z?^^-* 0. 
In the opposite limit ?/ — we obtain 



2l7r2 11 ^ / 2N 
+ — + [y^] 



32y 
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O a 



O a 



Thus, in the high-temperature regime t/k oo, we have 

,21n^r 



-Nra 



16 K 



(24) 
(25) 

(26) 



which indicates that this (3 function is becoming very large and negative. Similarly, as 
t/k — > CX) 



,21ti^t 



16 K 



(27) 



Our sign here is in agreement with that of Landsman 0. Unlike Landsman however we 
do not believe it to be an artefact of a one-loop approximation. The sign appears quite 
naturally as increasing temperature in a function that depends only on t/k is equivalent 
to decreasing k,, i.e. to going into the infrared. 

A calculation to two loops of the P function is expected to yield in the high-temperature 
limit a contribution proportional to a^/y"^, and higher- loop calculations would give higher 
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powers in temperature. In distinction to perturbation theory without RG improvement 
at the actual temperature T, increasing powers of 1/?/ in the high-temperature hmit are 
not necessarily a problem in our RG approach. This is seen by implementing the change 
of variables h = —agydf/dy |]lT|. One finds the expansion 

f3{h) = -e{y)h - + H{y)h^ + 0{h^), (28) 

where 

The function H{y) will be finite for < ?/ < oo. The function e{y) goes to 1 as ?/ — and 
vanishes in the limit ?/ — > oo. This expansion will be valid for all y as long as h is not too 
large. The transformation to /i is a change of dependent variable in a differential equation, 
it cannot change the physics. It does make manifest, however, the well behaved nature 
of the perturbative expansion of the (3 function for y ^ Q. In this limit the resultant (3 
function is exactly that expected of zero-temperature three-dimensional QCD. As y ^ oo 
one recovers the j3 function for zero-temperature four-dimensional QCD. In conclusion, as 
long as /i -C 1, the (5 function expansion, in particular its sign, is perturbatively reliable. 



5 Discussion 

For the moment we consider the number of fermions small enough to ensure asymptotic 
freedom at zero temperature, i.e. ^-Nc — \Nf > 0. At fixed, finite temperature, the 



running coupling ( ]22|) then vanishes in the high-momentum limit in much the same way 
as at zero temperature, as an inverse logarithm (see fig. 1 where iVc = 3 and Nf = 
6). The solutions of the differential equation show a pole at low momenta, where the 
coupling increases beyond the range of reliability of perturbation theory. Its position is 
governed by a rational function of the initial condition when k/T is small instead of a 
logarithmic one as in the zero-temperature case. The breakdown of perturbation theory 
in the infrared is not totally unexpected as there is no IR cutoff. No finite-temperature 
tree level magnetic mass was assumed and no magnetic mass is induced in one loop (even 
in RG improved perturbation theory as is used here). We note that for the running 
coupling, the low-momentum limit is not necessarily the same as the high-temperature 
limit, because changing the temperature may result in hopping to another curve with a 
different initial condition. This behaviour is determined by the temperature beta function 

The temperature-dependence of the coupling at fixed momentum scale is shown in 
fig. 2. At zero temperature the coupling approaches a finite value which can be taken 
to be the zero-temperature coupling at that scale. At the high-temperature end the pole 
shows up again. The existence of the pole can be seen by using the asymptotic expression 
(p7|). However to find its position in terms of zero-temperature parameters one needs 
the complete arbitrary-temperature expressions. The increase with temperature of the 
running coupling constant seems to contradict the notion of asymptotic freedom. However, 
the standard argument relies on a replacement of momentum scale by temperature 



in a vacuum result. This corresponds to the high-momentum limit discussed above. In 
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contrast, the coupling constant for growing T at fixed momentum scale increases. This 
is consistent with dimensional reduction |[13| which implies that at high temperature one 
is driven into the three-dimensional strong-coupling infrared regime. In view of ( ]28|) 
dimensional reduction takes place as ?/ — 0, where / /asymp := Nc21tt'^ /32y. However, 
this limit can be taken in many different ways. Also, in order for the result to have any 
validity the coupling constant in this limit must not be too large, say < 1. From 
the contour plots in fig. 3 one sees that as T/k increases, k itself must also increase in 
order that stays within the perturbative regime. More specifically, when / ^ /asymp, 
one requires that T/k increase no faster than {16/21n^){fNc - iiV/) ln(fi:/AQCD). Note 
that our contour plot is different from the one found in ref. |14] which shows a decreasing 
coupling. 

With the quark-gluon plasma in mind one particular trajectory suggests itself because 
there is, of course, a characteristic momentum scale for the system given by, for instance, 
the thermal average < p > of the momentum. For non-interacting, massless quanta 
in thermal equilibrium < p >= 3T. Following we may investigate the behaviour 
of the coupling constant along this trajectory. One finds the coupling to decrease as 
1/lnT. However, one should remember that this coupling is small only for fluctuations 
in the system which have a certain (albeit typical) momentum. For T decreasing one 
enters a strong coupling regime {as{n = 3T, T) > 1) at a temperature which depends 
on Nc and Nf. For A^^c = 3 and Nf = 6 perturbation theory breaks down for T less 
than Tb ~ I.SAqcd- Clearly in this region the original assumption of non- interacting 
particles is invalid. One might also consider as in |]T5[ an ensemble average of the coupling 
as{T) =< as{K,T) > by integrating over all momenta weighted by a statistical factor. 
However, then the problem is that the contributions to the average coupling from small 
momenta are too large to be trusted in the confines of a one-loop calculation, see fig. 
3. This point can be taken further to physical quantities in the magnetic sector of the 
equilibrium system. Such a quantity will be computed from an ensemble average over 
all fluctuation scales and it will be in the form of an expansion in as{T,p). Hence the 
contribution to the ensemble average from low-momentum fluctuations will be large and 
perturbatively unreliable. 

In the above we assumed that Nf < y-^c, in other words there are not enough fermions 
to prevent the colour charge from being screened at high energies. If, alternatively, Nf > 
Y^c the theory is no longer asymptotically free as k — > oo for fixed T. In fact the theory, 
at least perturbatively, is now ill-defined in that limit but is expected to be trivial. The 
interesting physics arises by noticing that as T/ k, becomes large, the fermions, because of 
the form of the fermionic distribution function, preferentially decouple from loops relative 
to the bosons. This can be seen from the expansions (|^) and (|25|) , the former being 
dominant relative to the latter in the y ^ limit. Although the fermions decouple 
from loops, this does not mean that the high-temperature behaviour is independent of 
the number of fermions. The coupling constant as(/«,T) is highly sensitive to the zero 
temperature initial conditions, i.e. as{K,, 0) which obviously depends on Nf. Observe that 
Jias a zero, which we take to he at y = yi (see fig. 4). In fig. 5 a contour plot 
is shown. At fixed, finite temperature, as k increases, decreases for k < 2Tyi and 
increases for k, > 2Tyi. On the other hand, for fixed k and along the line k = 3T the 
coupling monotonically increases as a function of T. 

So what are we to conclude from this one-loop calculation of a gauge-invariant coupling 
constant appropriate for the magnetic sector of QCD? One might expect that as in the 
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case of \(f)^ and QED one could completely control the high-temperature behaviour. For 
instance A0"^ is "strongly coupled" in the infrared at high temperatures but this we can 
deal with using RG techniques, since it has a non-trivial three-dimensional IR fixed point. 
For QCD however, we have found that for fixed momenta this couphng has a pole, thereby 
invalidating a RG-improved one-loop approximation. The reason for this is, we believe, 
that QCD is a confining theory — we simply do not know how the QCD coupling behaves 
in the confining regime. The strong coupling then is a result of dimensional reduction 
into the IR confining region of three dimensional QCD. Without being able to solve the 
confinement problem we obviously cannot track the crossover between quark-gluon degrees 
of freedom and confining degrees of freedom. It is quite natural to enquire as to whether 
this is a one loop artifact, especially in light of the fact that a magnetic mass is induced 
at two loops. Obviously this question cannot be answered without addressing higher loop 
effects. 
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Figure Captions 

Fig. 1: Plot of solutions of the P^'^^ differential equation versus k at fixed temperature 
for two different values of the integration constant. A^^c = 3 and Nf = 6. 

Fig. 2: Plot of solutions of the /^'-'^^ differential equation (|^ versus T at fixed momentum 
K for two different values of the integration constant. Nc = 3 and Nf = 6. 

Fig. 3: Contour plot of the running coupling constant ^^(/t, r) ( p2D for r = T and Tq = 
where Nc = 3 and Nf = 6. Only below the curve = oo (close to the curve = 1000) 
the coupling constant is positive and finite. Temperature and momentum are measured 
in units of Aqcd- The dashed line represents k = 3T. 

Fig. 4: Graph of the beta function P^'^\y) for different fermion numbers. For Nf > 
the beta function has a zero. 

Fig. 5: Contour plot of the running coupling constant as{n, r) ( P^D for r = T and Tq = 
where A^^c = 3 and Nf = 24. Temperature and momentum are measured in units of Aqcd- 



